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The most massive globular cluster in the Milky Way, u Centauri, is thought to be the re- 
maining core of a disrupted dwarf galaxy 1 ' 2 , as expected within the model of hierarchical 
merging 3 ' 4 . It contains several stellar populations having different heavy elemental abun- 
dances supplied by supernovae 5 — a process known as metal enrichment. Although M22 
appears to be similar to u Cen 6 , other peculiar globular clusters do not 7 ' 8 . Therefore u Cen 
and M22 are viewed as exceptional, and the presence of chemical inhomogeneities in other 
clusters is seen as 'pollution' from the intermediate-mass asymptotic-giant-branch stars ex- 
pected in normal globular clusters 9 . Here we report Ca abundances for seven globular clus- 
ters and compare them to uj Cen. Calcium and other heavy elements can only be supplied 
through numerous supernovae explosions of massive stars in these stellar systems 10 , but the 
gravitational potentials of the present-day clusters cannot preserve most of the ejecta from 
such explosions 11 . We conclude that these globular clusters, like u Cen, are most probably 
the relics of more massive primeval dwarf galaxies that merged and disrupted to form the 
proto- Galaxy. 

The Sejong/ ARCSEC Ca uvby survey program was initiated in 2006 to investigate the ho- 
mogenous metallicity scale for globular clusters and to obtain the complete metallicity distribution 
function of the Galactic bulge using the hk index [= (Ca — b) — (b — y)] u . The Ca filter in the hk 
index measures ionized calcium H and K lines, which have been frequently used to calibrate metal- 
licity scale for globular clusters 13 ' 14 . The utility of the hk index is that it is known to be about three 
times more sensitive to metallicity than the mi index is for stars more metal-poor than the Sun and 
it has half the sensitivity of the mi index to interstellar reddening 12 . During the last three years, 
we have used more than 85 nights of CTIO 1.0-m telescope time for this project. The telescope 
was equipped with an STA 4k x 4k CCD camera, providing a plate scale of 0.289 arcsec/pixel and 
a field of view of 20 x 20 arcmin. All of our targets accompanied with standards were observed 
under the photometric weather conditions and most of targets were repeatedly visited between 
separate runs. The photometry of our targets and standards were analyzed using DAOPHOT II, 
ALLSTAR, and ALLFRAME 1516 . 



In the course of metallicity calibration of red giant branch (RGB) stars in GCs, we found 
that many GCs show split in the RGB in their hk versus V color-magnitude diagrams (Figs 1 and 
2). The prime examples are M22 and NGC1851. In particular, the double RGB sequence in M22 
is very intriguing. The differential reddening effect and the contamination from the off cluster 
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populations cannot explain the double RGB sequences in M22 (see Supplementary Information). 
It has been debated for decades whether this cluster is chemically inhomogeneous or not, but the 
recent high resolution spectroscopic study of 17 RGB stars in the cluster suggests that it contains 
chemically inhomogeneous subpopulations 6 . The bimodality in the mi index of M22 RGB stars 
was also known, but it has been argued that it is most likely due to the bimodal CN abundances, 
where CN absorption strengths strongly affect the m x index, not due to the bimodal distribution 
of heavy elements in the cluster 1719 . The star-to-star light elemental abundance (C, N, O, Na, Mg 
and Al) variations have been known for decades and they are now generally believed to be resulted 
from chemical pollutions by intermediate-mass asymptotic giant branch stars 9 or fast rotating mas- 
sive stars 20 . However, it should be emphasized that our hk measurements for RGB stars in M22, 
NGC1851 and other GCs show discrete or bimodal distributions in calcium abundance, which 
cannot be supplied by intermediate-mass asymptotic giant branch stars or fast rotating massive 
stars. 

As shown in Fig 3, the difference in calcium, silicon, titanium and iron abundances between 
the calcium weak (Ca-w hereafter) group with smaller hk index and the calcium strong (Ca-s here- 
after) group with larger hk index in M22 and NGC1851 suggests that they are indeed chemically 
distinct 21 24 . (It is not shown in the figure but europium also has a bimodal abundance distribution 
in M22, in the sense that the Ca-s group has a higher europium abundance.) As for the origin of 
chemical inhomogeneity in globular clusters, at least four viable chemical enrichment mechanisms 
have been proposed up to date. They are, in the order of time required to emerge; (i) fast rotat- 
ing massive stars, (ii) Type II supernovae, (iii) intermediate-mass asymptotic giant branch stars, 
and (iv) Type la supernovae. If the current understanding of supernovae explosions is correct, 
only Type la and II supernovae can supply the heavy elements such as calcium and iron 10 . To ex- 
plain the discrete calcium abundances seen in M22 and NGC1851, however, the contribution from 
Type la supernovae can be ruled out for two reasons. First, the longer timescale (> 1 - 2 Gyr) 
before the onset of Type la supernova explosions, which would produce detectable age spread be- 
tween two populations; and second, the enhanced ^-elemental abundances, indicative of absence 
of contributions from Type la supernovae 10 . Qualitatively, the differences in elemental abundances 
between the two stellar populations in M22 and NGC1851 can be naturally explained by invoking 
chemical enrichment by Type II supernovae, where a-elements (silicon, calcium, and titanium) 
and r-process element (europium) are dominantly produced. However, our results do not neces- 
sarily imply that Type II supernovae are solely responsible for the chemical enrichment in M22 
and NGC1851, since all four above-mentioned mechanisms may be involved. We emphasize that 
the crux of our results is the undeniable evidence for Type II supernovae contribution to chemical 
enrichment of some globular clusters, in sharp contrast to the widely accepted idea of chemical 
pollution only by intermediate-mass asymptotic giant branch or fast rotating massive stars, with 
which the chemical enhancement of the a- and r-process elements in the second generation of the 
stars cannot be easily explained. 

More than half of 37 globular clusters in our sample shows discrete or broad distributions 
of the hk index in their RGB sequences. In Fig 2, We show color-magnitude diagrams for some 
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of exemplary globular clusters in the order of hk widths of RGB sequences at V HB , the V magni- 
tude level at the horizontal branch: u Cen, M22, NGC1851, NGC2808, M4, M5, NGC6752 and 
NGC6397 (see also Supplementary Table 3 and Figs 6 - 13). NGC2808 is known to have multiple 
main-sequences but no multiple RGB sequences have been reported to date. Our new results show 
that NGC2808 shows at least two discrete RGB sequences with a large spread in calcium abun- 
dance. Similarly, M5 has very broad hk index in the RGB sequence and NGC6752 shows discrete 
RGB sequences. It is interesting to note that all the globular clusters with signs of multiple stellar 
populations have relatively extended horizontal branch, while the globular clusters with normal 
horizontal branch (e.g. NGC6397 in Fig 2 and Supplementary Fig 13) show no spread or split in 
RGB. This is consistent with the suggestion that the extended horizontal branch is a signal of the 
presence of multiple stellar populations in globular clusters 25 . 

The overwhelming problem of the chemical enrichment by Type II supernovae in globular 
clusters is that their ejecta are considered to be too energetic to be retained by less massive systems 
like typical Galactic globular clusters (< a few times 10 5 M©) 11 . Our results therefore suggest that 
M22, NGC1851 and other globular clusters with RGB split were much more massive in the past, 
unless the current understanding of supernovae explosions is in great error. Perhaps, these globular 
clusters were once nuclei of dwarf-galaxy-like fragments and then accreted and dissolved in the 
Milky Way, as is widely accepted for cu Cen 1 ' 2 26 . Recent calculations suggest that a massive (> a 
few times 10 6 M ) star cluster embedded in a proto-dwarf galaxy could accrete gas from its host 
dwarf galaxy which may cause the formation of the second generation stars, producing multiple 
stellar populations 27 . Note that this scenario is also suggesting that the globular clusters with 
multiple stellar populations would be the remaining cores of the proto-galactic building blocks. 
This idea is supported by the recent investigations of the extended horizontal branch globular 
clusters (i.e. globular clusters with signatures of multiple stellar populations), which has shown that 
extended horizontal branch globular clusters are clearly distinct from the normal globular clusters 
in orbital kinematics and mass 25 . Extensive photometric surveys for fainter stars in these globular 
clusters, as well as spectroscopic surveys for stars in double RGB sequences, would undoubtedly 
help to shed more light into the discovery reported here. 
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Figure 1: Color-magnitude diagrams for M22. a, V versus b — y; b, V versus hk. In b, note the 
distinct and discrete double RGB sequences in M22. This cannot be due to differential reddening 
effect across the cluster or the contamination from the off cluster field but, is due to the difference 
in calcium abundance, which was synthesized in supernovae, between the two RGB sequences. 
The number ratio between the Ca-w group with smaller hk index and the Ca-s group with larger 
hk index is about 70:30. Black arrows in each panel denote reddening vectors. 
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Figure 2: Color-magnitude diagrams for lu Cen, M22, NGC1851, NGC2808, M4, M5, 
NGC6752 and NGC6397. Note that, while the distributions of the RGB sequences in the b - y 
color are relatively narrow, those in the hk index are either discrete or broad. This is evidence for 
the multiple stellar populations with distinct calcium abundances. Among these globular clusters, 
NGC6397 appears to be the only normal globular cluster with simple population (i.e. coeval and 
monometallic). 
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Figure 3: Differences in chemical compositions between double RGB sequences in M22 and 
NGC1851. a, b, Black 'plus' signs denote stars in M22 with proper motion membership prob- 
abilities P > 90%; blue filled diamonds and red filled circles denote RGB stars studied with 
high-resolution spectroscopy in the Ca-w and the Ca-s groups, respectively 23 24 . The green solid 
line denotes the fiducial sequence of RGB stars and Ahk denotes the difference in the hk index 
against the fiducial sequence. The double RGB sequences persist in proper motion member stars, 
c - f, Comparisons of elemental abundances between the Ca-w and the Ca-s groups in M22. Solid 
lines denote the mean values, and dashed lines denote standard deviations of each group. The Ca-s 
group has higher a-elements (Si, Ca, and Ti) and iron abundances, which must be supplied by 
numerous Type II supernova explosions, g, h, Black 'plus' signs denote stars in NGC1851; blue 
filled diamonds and red filled circles denote RGB stars studied with high-resolution spectroscopy 
in the Ca-w and the Ca-s groups, respectively 2122 , i - 1, As c - f but for NGC1851. 
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Supplementary Information 



1 The CTIO Ca Filter System 

The Filter Transmission The Ca filter system was designed to include Ca II H and K lines at 
A 3968 and 3933 A, respectively, with a full-width half maximum (FWHM) of approximately 90 
A. The CN band absorption strengths at ps A 3885 A are often very strong in stellar spectra and 
the lower limit of the Ca filter is set to avoid contamination by the CN band 12 . The Ca filter used 
at Cerro Tololo Inter-American Observatory (CTIO) has a similar FWHM, approximately 90 A, 
but its passband is shifted approximately 15 A to the longer wavelength (Alistair Walker, private 
communication) compared to that in Anthony-Twarog etal. 12 In Supplementary Figured! we show 
the transmission functions of the both Ca filters. In the figure, we also show synthetic spectra for 
the CN normal and the CN strong RGB stars in typical intermediate metallicity globular clusters 
(GCs) as an illustration. The effect of the CN band on the hk index is negligible as will be discussed 
below. 

The Central Wavelength Drift of the Ca Filter The CTIO Ca filter transmission function shown 
in Supplementary Figure[T]is that measured with a collimated beam. It is known that the passband 
of the narrow band interference filter depends on the angle of the incidence beam following, 



where Ao is the wavelength of peak transmittance at normal incidence, (3 is the angle of incidence 
of the collimated beam on the filter and n* is the effective refractive index of the filter 28 . Therefore, 
when the Ca filter is used with a fast telescope, the filter passband can be significantly different 
from that shown in Supplementary FigureQ] The CTIO 1-m telescope used for our survey is a slow 
telescope with //10.5 and the effect resulted from the angular dependency of a converging beam is 
expected to be very small. Assuming [3 ~ 1/21 radians for the converging beam at the CTIO 1-m 
telescope and n* « 1 .4 for the CTIO Ca filter, the peak wavelength of the Ca filter will be shifted 
by 2.3 A to the shorter wavelength. Given the much larger FWHM of the CTIO Ca filter, this may 
contribute small effect. We investigate contributions to the hk index resulted from the shifted Ca 
passband using synthetic spectra for typical intermediate metallicity RGB stars in our GCs. Our 
calculations integrating over the filter transmission curve show that this effect contributes no larger 
than 0.011 mag to our hk measurements, in the sense that the shifted Ca passband to the shorter 
wavelength gives slightly larger hk values. We emphasize that, since our results are based on a 
single instrument setup (the same telescope, filters and the CCD camera) during the observations 
of our science targets and the photometric standards, this effect is expected to be cancelled out 
during our photometric calibrations. Also importantly, our main results presented here rely on the 
split or the spread in the hk index of RGB stars of an individual GC. Therefore, the shifted Ca 
passband affects similar degree to the hk index of the RGB stars in a GC and does not contribute 
to the apparent RGB split or the spread in the hk index of an individual GC. 




(1) 
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Effect of radial motions and internal velocity dispersions of GCs The mean radial motion of 
139 GCs in our Milky Way Galaxy 29 is |u r | = 1 10 km/s, equivalent to the wavelength shift by | AA| 
~ 1 .4 A at A 3950A. We calculate the contribution due to the mean radial motion of GCs to the hk 
index using the shifted CTIO Ca passband and synthetic spectra for typical intermediate metallicity 
GC RGB stars. We find that the net effect is negligibly small, \Ahk\ < 0.006 mag. Among our 
eight GCs, NGC1851 has the largest radial velocity, v r = 321 km/s, equivalent to the wavelength 
shift by AA = 4.2 A to the longer wavelength at A 3950A. We calculate the contribution due to 
the radial motion of NGC1851 using the shifted CTIO Ca passband and the red-shifted synthetic 
spectra with a fixed CN abundance for the cluster. We obtain Ahk ~ 0.015 mag, in the sense that 
the red-shift is resulted in a slightly larger hk index. As we discussed above, the difference in the 
hk index due to the high radial motion of NGC1851 does not affect our results presented here, 
since the hk indices of RGB stars in NGC1851 will be affected by similar degree and the mean 
radial motion of the cluster does not produce an apparent split or a spread in the hk index. Perhaps, 
this effect may become important in the inter-cluster comparisons, which is beyond the scope of 
our study. 

What concerns us most about the high radial velocities of some GCs, in particular for red- 
shift, is the potential contamination by the strong CN band absorption features at A 3885 A as 
shown in Supplementary Figured! For example, NGC1851 has a bimodal CN distribution and 
some RGB stars show very strong CN band absorption strengths 30 . Due to its high radial velocity 
away from us (i.e. red-shifted), the CN band absorption features in the CN-strong RGB stars may 
affect the hk index and, subsequently, may produce an apparent RGB split of the cluster as shown 
in Figure 2 or Supplementary Figure [8l We calculate the CN band contributions using the shifted 
CTIO Ca passband and the red-shifted synthetic spectra for the CN-normal and the CN-strong RGB 
stars (see discussion below). Our calculations suggest that the net effect is negligibly small, Ahk 
< 0.003 mag, and the high radial velocity of NGC 1851 combined with a bimodal CN distributions 
does not produce the RGB split in the hk index. 

We also investigate the effect of the internal velocity dispersion of an individual GC. As- 
suming o~los = 15 km/s, equivalent to AA < 0.2 A at A 3950A, we obtain Ahk < 0.001 mag 
following the same method described above, and the effect from the internal velocity dispersions 
of GCs does not affect our results. 

Summary of uncertainties on the hk index Supplementary Table \T\ summarizes the uncertain- 
ties in our hk index measurements relevant to the CTIO Ca passband. (The variations in the hk 
index due to differences in elemental abundances of GC RGB stars will be discussed below.) As 
discussed above, the effects due to the shifted CTIO Ca passband and the radial motions of GCs 
do not affect our results, since both effects contribute similar degree to the hk indices among RGB 
stars in a GC. (i.e. They only affect the zero point of the hk index and they do not affect the 
Ahk distributions). In addition to our photometric measurement errors which will be discussed 
below, the effects due to the internal velocity dispersions of GCs and the differential interstellar Ca 
II absorption (see discussion below) can affect our hk index measurements. However, their con- 
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tributions to our hk measurements are no larger than 0.022 mag and they do not affect our main 
conclusion presented here. Therefore, our results strongly suggest that the split or the spread in 
the hk index of RGB stars in GCs are related to the variations in elemental, in particular calcium, 
abundances among RGB stars in a GC, which will be discussed below. 

2 The Double RGB Sequences of M22 

Differential Reddening Effect on the Double RGB Sequences in M22 The continuous interstel- 
lar extinction by the interstellar dust and the discrete interstellar line extinction by the interstellar 
Ca II atoms may affect our main results. We considered both effects and will discuss that the RGB 
split of M22 in the hk index is indeed due to the difference in calcium abundances between two 
stellar populations in M22 and other explanations are highly unlikely. Also both effects tend to 
produce spreads in RGB sequences rather than the distinct and discrete RGB sequences of GCs 
reported here. 

The differential continuous reddening across the cluster can thicken the apparent RGB se- 
quence of GCs in broad-band optical photometry 31 . In contrast to other color indices being used 
in broad-band photometry, the hk index is known to be insensitive to interstellar reddening 12 , 
E(hk) /E(b -y) = -0.16 and E(hk) /E(B -V) = -0.12. The difference in the hk index between 
the two RGB sequences in M22 is about 0.2 mag at the magnitude level of the horizontal branch. 
If this hk split is only due to differential reddening effect, we would expect even larger separation 
of the two RGB sequences in the b — y color and the V magnitude. The reddening correction value 
in the b — y color, E(b — y), for the Ca-s group is about —1.25 mag, equivalent to E[B — V) = 
— 1.69 mag assuming E(b — y)/E(B — V) = 0.74, making the RGB stars in the Ca-s group too 
hot to be RGB stars (see Supplementary Figure [2]- d). At the same time, the extinction correction 
value in the V magnitude is —5.24 mag, assuming Ay = 3AxE(B — V), for the Ca-s RGB stars. 
Applying this large extinction correction makes the RGB stars in the Ca-s group too bright to be 
members of M22 (see Supplementary Figure[2]- e & f). We emphasize that both the Ca-w and the 
Ca-s groups are proper motion members of the cluster as shown in Figure 3. Note also that the 
reddening vector (see Figure 1 or Supplementary Figure 0) is almost parallel to the slopes of HB 
and RGB in the hk versus V CMD, and thus the differential reddening can not produce the RGB 
split. Therefore, continuous differential reddening effect can be completely ruled out to explain the 
observed bimodal RGB sequences in M22. Similarly, the interstellar reddening toward NGC1851 
is very small, E(B — V) = 0.02 mag 29 , but the hk split in RGB stars of the cluster is as large 
as 0.2 mag (see Supplementary Figure [8]), which can not be explained by differential continuous 
reddening effect 22 . 

The previous study for the GCs showed that the equivalent width of the interstellar Ca II K 
absorption line strength can be as large as several times 100 mA 32 . The interstellar Ca II atom is 
thought to be heavily depleted on to dust in denser clouds 33 , which may cause small-scale differ- 
ential discrete reddening effect across M22 and other GCs studied here (see also Andrew et al. u 
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for the small-scale variations of interstellar Na I D lines^] toward the less extincted globular cluster 
M92). We generate synthetic spectrum to surrogate interstellar Ca II H & K absorption lines. We 
adopt a gaussian line profile with a FWHM of 1 A, equivalent to Av r ~ 76 km/s, and we assign 
equivalent widths of 350 mA and 650 mA for the interstellar Ca II H & K lines, respectively. Our 
synthetic spectrum is shown in Supplementary FigureQ]- (c). Assuming they are linear part of the 
curve of growth 35 , the column density of the interstellar Ca II can be estimated as 

JV(Call) = 1.13 x 10 20 — -, (2) 

A 2 / 

where EW and A are the equivalent width and wavelength in A and / is the oscillator strength. 
Using the oscillator strengths of 0.681 and 0.341 for Ca II H & K, respectively, the column density 
for interstellar Ca II is log iV(Ca II) « 12.8 cm" 2 , equivalent to AE(B - V) » 0.32 mag 33 . If this 
large amount of small-scale interstellar Ca II variation exists among our GCs studied here, how 
much will it affect our hk index measurements? We calculate the Ca magnitudes with and without 
the interstellar Ca II variations using the shifted CTIO Ca transmission function. The difference in 
the Ca magnitude (i.e. in the hk index since the interstellar Ca II H & K lines do not affect b or 
y passbands) is only 0.010 mag and, therefore, the differential discrete reddening effect due to the 
variations in the interstellar Ca II abundances can be completely ruled out to explain the GC RGB 
splits in the hk index. 

The Spatial Distributions In Supplementary Figure [31 we show the spatial distributions of Ca-w 
and Ca-s RGB stars in M22. As can be seen in the figure, each population does not show any 
spatially patched features, supporting our results that differential reddening is not responsible for 
the RGB split in M22. 

Contamination from the Milky Way's Bulge Population M22 is located in the direction of 
the Milky Way's bulge and the contamination from the bulge population may affect our results. 
However, this is very unlikely, since the proper motion member RGB stars show discrete double 
RGB sequences as shown in Figure 3. In addition, the bulge RGB stars are located farther from the 
Sun, more metal-rich and suffering from heavier interstellar reddening than those in M22 are. In 
Supplementary Figure HI we compare M22 CMDs with those of two bulge fields (NGC6528 and 
OGLEII - 12). As can be seen in the figure, the RGB stars in the bulge are fainter and redder than 
those in M22 are and the contamination from the Milky Way's bulge population does not affect 
our results. 

Effects of Metal Contents and Helium Abundances on the M22 RGB As shown in Figure 3, 
the stars in the Ca-s group are about 0.2 dex more metal-rich than those in the Ca-w group. It 
is suspected that this large metallicity spread may produce any detectable discrepancy in stellar 
evolutionary sequences, in particular for RGB sequence, between two stellar populations based 

'Note that the number of interstellar Na I atoms appears to be about a factor of ten larger than that of interstellar 
Ca II atoms 33 . 
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on broad-band photometry. To explore metallicity effect on the RGB sequence, we compare BV 
CMD by Monaco et al? 1 with the latest Y 2 isochrones (Version 3, Yi et al. in preparation). In 
Supplementary Figure [5J we show model isochrones for [Fe/H] = — 1.6 and —1.4 with the helium 
abundance of Y = 0.23 and the age of 1 1 Gyr, using the reddening value and the distance modulus 
for the cluster from Harris 29 . Although the split in the RGB sequences of two model isochrones 
is noticeable, the discrepancy in the RGB sequence does not appear to cause a serious problem to 
explain the BV CMD by Monaco et al? 1 Note that the V magnitude difference in the sub-giant 
branch between two model isochrones can be as large as 0.2 mag, apparently consistent with recent 
HST/ACS observations of the cluster 36 . 

As inferred from the extended HB (EHB) morphology of M22, the second generation of the 
stars is expected to have enhanced helium abundance by AY « 0.05 37 . Since the new version 
of Y 2 isochrones provides models with enhanced helium abundances, we investigate the effect 
of helium abundance on the evolutionary sequence. As illustrated in Supplementary Figure [5]- 
(c), the discrepancy between two stellar populations alleviates due to the opposite effect of metal 
contents and helium abundances on the RGB temperature. Since the second generation of the stars 
in M22 shows signs of the chemical enrichment by Type II supernovae and intermediate-mass 
asymptotic giant branch (AGB) stars, the second generation of stars may be slightly younger than 
the first generation. Assuming the age difference of 1 Gyr between two stellar populations, two 
model isochrones are in excellent agreement except for bright RGB sequence, where the observed 
number of stars is small. 

It is intriguing to note that the number ratio between the Ca-w and the Ca-s RGB stars (70:30) 
found here is very similar to those found between (1) the two stellar groups with different [Fe/H] 
and [s-process/Fe] ratios 6 , (2) the brighter SGB and the fainter SGB stars 36 , and (3) the two groups 
of HB stars with the bluer HB being less populated. The population synthesis models (Han et 
al. 2009, in preparation) suggest that this can be naturally reproduced by the enhanced metal and 
helium abundances in the second generation of stars. 

3 The hk and Metallicity Distributions of GCs 

Observations In Supplementary Tabled we show the journal of observations for eight GCs. They 
were observed under the photometric weather conditions and, for most cases, the median seeing 
was about 1.5 - 1.6 arcsec during our observations. Note that the RGB stars in NGC2808 are 
roughly 21 times fainter than those in NGC6397 in the Ca passband for a fixed magnitude, while 
our total integration time for NGC2808 is only about three times longer than that for NGC6397 in 
the Ca passband. Statistically, the lack of total integration time for NGC2808 will be resulted in a ~ 
2.6 times larger Ca measurement error than that expected for NGC6397 at a given Ca magnitude, 
particularly for fainter stars. Although our survey relied on a rather small telescope through a rather 
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narrow filter at a rather blue wavelength^, our investigations of the multiple stellar populations of 
GCs are focused on bright RGB stars, where the numbers of photon in the Ca passband are enough 
so that the measurement errors, including propagation errors during the photometric calibrations, 
are less than 0.020 mag (see Supplementary Table [3]). 

Color Distributions of Bright RGB Stars Here, we investigate the b — y color and the hk index 
distributions of RGB stars brighter than V — Vhb = 1.0 mag. We derive lower order 4 - 5) 
polynomial fits, which are forced to pass through the peak b — y colors or the peak hk indices of 
given magnitude bins, as fiducial sequences for eight GCs and then we calculate differences in the 
b — y color, A (b — y), or the hk index, Ahk, with respect to fiducial sequences of each GC. We 
show our results in Supplementary Figures [6]— |T3} In the figures, the blue horizontal bars denote the 
mean measurement errors including propagation errors during the photometric calibrations with a 
2cr* range (± ler*) for individual stars at given magnitude bins. 

In Supplementary Table [H we show comparisons of the observed FWHMs of RGB stars 
and the measurement errors at the magnitude level of horizontal branch stars, Vhb + 0.5 > V > 
Vhb — 0.5, in each GC. To calculate the FWHMs of RGB stars in GCs, we used the following 
relation, 

FWHM(RGB) « 2.3548 x a A , (3) 

where is the standard deviation of RGB stars in the A(b — y) or the Ahk distributions in a 
fixed magnitude bin. Note that the FWHM of RGB stars in the hk index is slightly larger than 
the hk separation between the double populations such as in M22, NGC1851, NGC2808, M4 and 
NGC6752. The mean measurement errors, a*(b — y) or a*(hk), given in the table include propa- 
gation errors during the photometric calibrations and they are those for individual stars. Therefore, 
the mean measurement errors for an individual population in GCs, a p (b — y) or a p (hk), will be 
given by « a*(b — y)I^Jn^ or a*(hk)/ \frip~, where n p (> 20) is the number of stars in each popu- 
lation. Note that, while the FWHMs of most GCs have much larger values [> 8cr*(/iA;)] than the 
measurement errors for individual stars in the hk index, the FWHM of NGC6397 is comparable in 
size to the measurement error in the hk index, consistent with the idea that NGC6397 is the only 
normal GC with a simple stellar population (i.e. coeval and monometallic) in our sample. Also 
note that NGC6397 shows similar degree of the full RGB widths in the A(b — y) and the Ahk 
distributions (see Supplementary Tableland Supplementary Figure [TSl. 

The last two columns of Supplementary Tabled E{b — y)i/2 and E(hk)i/2, are for the 
contributions due to the continuous differential reddening effect assuming a 50% variation in the 
total interstellar reddening across each GC. The observed FWHMs of GCs in the (b — y) color 
and in the hk index can not be explained simultaneously, similar to what shown in Supplementray 
Figure [2l Therefore, the continuous differential reddening effect can be ruled out to explain the 

2 Fortunately, the CCD camera used in our survey has rather high quantum efficiency (QE) at shorter 
wavelength with QE w 0.686 at A 3800A and w 0.770 at A 4000A. (see http://www.astronomy.ohio- 
state.edu/Y4KCam/OSU4K/index.html#DQE). 
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differences between the observed FWHMs(RGB) and the measurement errors in GCs. 

Ahk as a probe of multiple stellar populations in GCs For u> Cen, M22 and NGC1851 (see 
Figure 3 and Supplementary Figure [T5l), when the two subpopulations are defined by our hk index 
(or our Ahk distribution), we can also see the clear division in spectroscopic elemental abundances. 
Similarly, we will discuss that the split or the spread in the Ahk distributions of RGB stars in other 
clusters can provide a powerful method to probe the multiple stellar populations in GCs. 

The calcium abundance is the major factor that determines the hk index or the Ahk distri- 
bution of RGB stars in a GC (see discussion below) and Type II supernovae are responsible for 
the calcium enrichment in a GC. As discussed, however, our results do not imply that Type II 
supernovae are solely responsible for the chemical enrichment in GCs. In an attempt to explain 
the observed large star-to-star lighter elemental abundance variations (in particular O and Na) in 
GCs, chemical pollution by intermediate mass AGB stars 9 or fast rotating massive (FRM) stars 20 
has been widely accepted. It should be reminded that, however, neither AGB nor FRM scenarios 
can explain the chemical enrichment of the a- and r-process elements in the second generation of 
the stars. It is most likely that all three aforementioned mechanisms (and perhaps including Type 
la supernovae) are required to explain elemental abundance patterns found in GCs. In addition 
to the chemical enrichment by Type II supernovae, which is the main results presented here, if 
the second generation of the stars in some of our GCs have experienced the chemical pollution 
by intermediate-mass AGB or FRM stars, the lighter elemental abundances, such as oxygen and 
sodium, between the two generatrions of stars must have been different. Furthermore, the vari- 
ations in [O/Fe] and [Na/Fe] can be as large as 1 dex in some GCs 7 and the differences in the 
oxygen and sodium abundances are easily detectable compared to those in the heavy elements, 
such as calcium and iron. 

During the last few years, tremendous amount of effort has been directed at spectroscopic 
study of RGB stars in GCs, in particular, using the multi-object spectrograph mounted at VLT. 
Among our eight GCs, NGC2808 38 , M4 39 and NGC6752 40 have been studied using moderately 
high resolution spectra for more than 100 RGB stars. In Supplementary Figure[14l we show com- 
parisons of Ahk versus O, Na and Fe distributions of the clusters. In panel (a), we show the plot 
of V - Vhb versus Ahk for NGC2808 RGB stars. In the figure, the plus signs denote the RGB 
stars with known [O/Fe] and [Na/Fe] ratios 38 . From the Ahk distribution of RGB stars shown 
in panel (b), we define the boundary at Ahk = —0.05 mag (the vertical dashed line) assuming 
that NGC2808 has two major stellar populations as shown in panel (b) or Supplementary Figure|9l 
Similar to the procedure employed in M22 and NGC 1851 (see Figure 3), we define the Ca-w group 
with smaller hk index and the Ca-s group with larger hk index and they are denoted by the blue 
and the red plus signs, respectively, in panel (a). In panels (c), (d) and (e), we show the [O/Fe], 
[Na/Fe] and [Fe/H] distributions for each group, where the shaded histograms outlined with blue 
color are for the Ca-w group and the blank histrograms outlined with red color are for the Ca-s 
group. The Ca-w group has a higher mean oxygen and a lower mean sodium abundances, while 
the Ca-s group has a lower mean oxygen and a higher mean sodium abundances, indicative of 
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the presence of the proton-capture process at high temperature between the two formation epochs 
presumably via intermediate-mass AGB or FRM stars, where oxygen is depleted by the CNO cy- 
cle while sodium is enriched from the 22 Ne + *H — > 23 Na reaction. Our results strongly suggest 
that they are truly different stellar populations and not related to, for example, our photometric 
measurement errors and differential reddening effect: the Ca-w group is the first generation of 
stars while the Ca-s group is the second generation of stars enriched by Type II supernovae (e.g. 
calcium) and intermediate-mass AGB or FRM stars (e.g. sodium). Although the difference in the 
[Fe/H] distributions between the two groups does not appear to be as compelling as those in the 
[O/Fe] and the [Na/Fe] distributions, the Ca-w group has a slightly lower mean metallicity than 
the Ca-s group does. We performed non-parametric Kolmogorov-Smirnov (K-S) tests to see if the 
[Fe/H] distributions of the two populations in NGC2808 are drawn from the same parent popula- 
tion. Our calculation shows that the probability of being drawn from identical stellar populations 
is 5.5% for NGC2808, suggesting that they have different parent populations. 

The same results can be found in M4 and NGC6752. From the comparisons of the [O/Fe] 
and the [Na/Fe] distributions between the Ca-w and the Ca-s groups, it can be seen that the Ca-w 
groups are the first generations of stars while the Ca-s groups are the second generations of stars in 
the clusters. We also performed K-S tests for the [Fe/H] distribution of M4, we obtained that the 
probability of being drawn from identical parent populations is 5.5% for M4, indicating that each 
subpopulation in M4 has different parent populations. 

Recalibration of [Fe/H]^ Based on RGB Stars in u Cen and Metallicity Distributions of 
Eight GCs In our previous study for NGC1851, we showed that the hk index traces the calcium 
abundance and, furthermore, it can provide a very powerful method to distinguish multiple stellar 
populations in GCs 22 . However, it can be seen that the full range of Ahk increases with the 
luminosity of RGB stars (i.e. different temperature or surface gravity), in particular, in to Cen 
and M22. Due to the temperature dependency on the hk index versus metallicity relation, the 
Ahk distributions cannot be directly translated into the absolute metallicity scale. Therefore, we 
calculate the photometric metallicity, [Fe/H]/^., of individual RGB stars in eight GCs using the 
[Fe/H] relations on the hko versus (b — y)o plane 12,22 . 

Recently, Johnson et al} studied elemental abundances, including calcium, of large sam- 
ple of RGB stars in uj Cen using moderately high resolution spectra (R ~ 18,000). Since uj Cen 
contains multiple stellar populations with very broad metallicity range, A [Fe/H] 1.5 dex, com- 
parisons of our results of RGB stars in uj Cen with those of Johnson et al. may provide an won- 
derful opportunity to assess our photometric metallicity scale using the hk index, [Fe/H] ftfc . In 
Supplementary Figure [T51 we show elemental abundances of 40 RGB stars in u Cen studied by 
Johnson et al. as a function of Ahk. As shown in the figure, [Ca/H] and [Fe/H] appear to be well 
correlated with Ahk, indicating that Ahk can truly be treated as the relative calcium abundance 
or metallicity indicators for RGB stars with similar luminosities in a GC. We also show plots of 
[Fe/H] spcc versus [Fe/H]/^ and [Ca/H] spcc versus [Fe/H]/^ for 32 RGB stars with sufficiently high 
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signal-to-noise ratios (> 100). We derive linear fits to each relation and we find 



[Fe/H] 



spec 



0.533 [Fe/H] 



hk 



0.775 (a = 0.087dex), 



(4) 



and 



[Ca/H] 



spec 



0.587[Fe/H] 



hk 



0.403 (a = 0.106dex). 



(5) 



We recalibrate our photometric metallicity using the equation ([4]), [Fe/H]/ l fe corr , and we de- 
rive metallicity distribution functions (MDFs) for eight GCs. During our calculations of MDFs, 
we use RGB stars with —2.0 < V — Vhb < —0.5 mag in order to minimize contamination from 
off-cluster field and red-clump populations. We show our results in Supplementary Figure PT5l As 
expected from the Ahk distributions, the signs of multiple stellar populations persist in our MDFs 
for most GCs. 

Finally, cautions are advisable on our MDFs of GCs. Our metallicity scale is not on the 
traditional Zinn & West 14 scale, therefore our MDFs for GCs can be different from those from 
other photometric or spectroscopic studies. RGB stars in uj Cen of Johnson et al. have different 
individual elemental abundances, which were not taken into consideration in our [Fe/H] spcc versus 
[Fe/H]h/c or [Ca/H] spec versus [Fe/H]/^. relations. Furthermore, each GCs may have slightly dif- 
ferent elemental abundance ratios and our calibrated photometric metallicities for GCs would be 
affected. However, it should be emphasized that the crux of our results is the split or the spread 
in the hk index in the RGB stars of individual GCs, which is insensitive to other elemental abun- 
dances except calcium as will be discussed below. 

4 The Influence of Elemental Abundances on the hk Index 

The realistic modeling of the resonance Ca II H & K lines requires proper understanding of stellar 
atmospheres, including chromospheres, and hydrodynamic non-local thermodynamic equilibrium 
treatments, which have posed difficult problems for decades 41 . Here, we demonstrate that the cal- 
cium abundance is the major factor that determines the hk index of RGB stars using 1 -dimensional 
plane-parallel stellar atmospheres 42 . 

Calcium Using the model atmosphere for the RGB star at the magnitude level of the horizontal 
branch with T cS = 4750 K, log g = 2.0 (in cgs unit), v turh = 2.0 km/s, [Fe/H] = —1.6, we calculate 
synthetic spectra for [Ca/Fe] = 0.25, 0.30, 0.35, 0.40, 0.45 and we show some of our synthetic 
spectra in Supplementary Figure [FT] We convolve the filter transmission functions with synthetic 
spectra and we obtain the calcium abundance sensitivity on the hk index, d(hk) /<9[Ca/H] « 0.422 
mag/dex. Note that this result is based on the fixed model parameters, such as T e g, log g, t>t ur b> an d 
[Fe/H], except calcium abundance. To interpret observed Ahk between the two stellar populations 
in a GC in terms of different calcium abundances, proper atmospheric parameters should be taken 
into consideration. For example, the two stellar groups with different [Fe/H] and [s-process/Fe] 
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ratios in M22 by Marino et al. have slightly different elemental abundances and temperatures. 
The stars in the metal-poor group by Marino et al. have ([Fe/H]) = —1.82, ([Ca/Fe]) = +0.25 and 
those in the metal-rich group have ([Fe/H]) = —1.68, ([Ca/Fe]) = +0.35. The mean temperature of 
the stars in the metal-poor group is ~ 100 ± 42 K hotter than those in the metal-rich group. Using 
these atmospheric parameters, we obtain Ahk = 0.121 ± 0.072 mag, which is apparently consistent 
with the double peaks in the Ahk distribution of M22 within the error as shown in Supplementary 
Figure [7] 

Helium As shown in Supplementary Figure \5\ helium is very important in stellar structure and 
evolution. Very unfortunately, however, there is no direct method to measure helium abundances of 
stars in GCs. As we discussed, all the GCs with signs of multiple stellar populations have relatively 
EHB, for example, the second generation of the stars in M22 is expected to have enhanced helium 
abundance by AY « 0.05 inferred from its EHB morphology. Using the model atmospheres with 
enhanced helium abundance (Y m 0.35, equivalent to AY « 0.10) by Castellf], we obtain Ahk « 
—0.002 mag, in the sense that the hk index decreases as helium abundance increases, and thus the 
effect of enhanced helium abundance on the hk index is negligible. Given the cool temperatures of 
RGB stars in GCs, the helium enhancement by AY = 0.05 - 0.10 does not appear to be important. 

CNO It is well-known fact that many GCs show large star-to-star elemental abundance variations. 
In particular, almost all GCs show variations in the CNO abundances resulted from the internal 
evolutionary mixing accompanied with the CNO-cycle or the primordial pollution by intermediate- 
mass AGB stars to the second generation of the stars 43 ' 44 . In spite of their high abundances, the 
CNO abundances are hard to measure in the optical wavelength mainly due to the lack of atomic 
transitions. On the other hand, in the form of molecules, the CNO can affect the hk index, in 
particular the CN band at A 3885 A as shown in Supplementary Figure [Q Both carbon and 
nitrogen contribute in the formation of CN molecules. The typical RGB stars in GCs show an 
anticorrelation between the CN band and the CH band strengths and a correlation between the CN 
band and the NH band strengths, indicating that the nitrogen controls the CN band strength 45 . Our 
results suggest that the variations in the CNO abundances do not affect the hk index significantly. 
We obtain Ahk « -0.007, +0.002, -0.004 for A[C,N,0/Fe] = +1.0 dex, respectively, and their 
influence on the hk index appears to be negligible. 

Aluminium It is also well-known fact that many GCs show large star-to-star aluminium variations 
by more than A[A1/Fe] ^1.0 dex, presumably resulted from the proton-capture process at high 
temperature or the primordial pollution by intermediate-mass AGB stars to the second generation 
of the stars 44 ' 46 . The resonance lines of Al I at A 3944.01 and 3961.52 A are often very strong 
(see Supplementary Figure ITTb and it may affect our conclusions that the hk index traces calcium 
abundances of RGB stars in GCs. We obtain Ahk w 0.013 mag for A[A1/Fe] = +1.0 dex. The 
effect of the variations in aluminium abundances on the hk index is insignificant compared to our 
observations, by more than a factor of ten. The insignificant influence of aluminium on the Ca II 

3 http://wwwuser.oat. ts.astro.it/castelli/grids. html 
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H & K lines was also confirmed by others for the globular cluster NGC6752 7 . 

a-elements The a-elements (O, Ne, Mg, Si, S, Ar, Ca, and Ti) are quite abundant and they are 
major donors to the H opacity in RGB stars. However, their influence, except for Ca, on the hk 
index appears to be small. We obtained Ahk fa —0.008 mag for +0.3 dex variation in a-elements 
excluding calcium. 

s-process elements RGB stars in GCs show large star-to-star s-process elemental abundance vari- 
ations presumably resulted from the primordial pollution by intermediate-mass AGB stars to the 
second generation of the stars 9 ' 44 . For example, globular clusters M22 and NGC1851 show bi- 
modal s-process elemental abundance distributions with A[s-process/Fe] > 0.5 dex. We obtained 
Ahk « +0.008 mag for A[s-process/Fe] = +0.5 dex, and thus their influence on the hk index is 
small. 
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Uncertainty on the hk index 


Note 


Photometry 


< 0.020 mag 


random 


Shifted CTIO Ca passband 


< 0.011 mag 


systematic 


Mean radial motion of GCs 


< 0.006 mag 


systematic 


Internal velocity dispersion 


< 0.001 mag 


random 


Interstellar Ca II absorption 


< 0.010 mag 


random 


total 


< 0.024 mag 




total (random) 


< 0.022 mag 





Supplementary Table 1 : Summary of uncertainties relevant to the hk index. 
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ID V HB E(B-V) Exposure Time (sec) Obs. Pos. Date 

Ca u v b y RA DEC (MIWYY) 



u> Cen 


14.53 


0.12 


12,860 




4,890 


2,130 


1,300 


13:26:44 


-47:26:28 


05/07, 02/08 


M22 


14.15 


0.34 


8,100 


2,400 


1,200 


2,530 


1,500 


18:36:29 


-23:55:34 


07/08, 08/08 


NGC1851 


16.09 


0.02 


19,100 


12,300 


7,400 


7,100 


3,810 


5:14:14 


-40:01 :49 


02/08, 08/08 


NGC2808 


16.22 


0.22 


10,820 




3,600 


4,960 


3,080 


9:11:57 


-64:49:24 


05/07 


M4 


13.45 


0.36 


8,400 


5,400 


5,570 


2,920 


2,060 


16:23:33 


-26:30:47 


05/07, 08/08 


M5 


15.07 


0.03 


9,660 


3,900 


2,100 


4,010 


2,390 


15:18:29 


2:04:03 


05/07, 08/08 


NGC6752 


13.70 


0.04 


7,500 


2,400 


1,800 


2,400 


1,200 


19:10:57 


-60:00:20 


07/08, 08/08 


NGC6397 


12.87 


0.18 


3,560 


3,560 


2,140 


1,355 


930 


17:40:52 


-53:36:06 


08/06, 09/07 



Supplementary Table 2: Journal of observations for eight GCs. Only one field has been 
observed for a particular GC and the coordinates are given in columns (9) and (10). 
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ID 


FWHM(RGB) 


Measurement errors 


Differential Reddening 


(b-y) 


hk 


o* (b - y) 




E(b-y)i/ 2 


E(hk) 1/2 


u) Cen 


0.079 


0.534 


0.012 


0.020 


0.089 


0.014 


M22 


0.050 


0.216 


0.004 


0.008 


0.252 


0.041 


NGC1851 


0.035 


0.182 


0.006 


0.013 


0.015 


0.002 


NGC2808 


0.042 


0.159 


0.008 


0.019 


0.163 


0.026 


M4 


0.037 


0.119 


0.005 


0.010 


0.266 


0.043 


M5 


0.025 


0.105 


0.006 


0.013 


0.022 


0.004 


NGC6752 


0.022 


0.090 


0.004 


0.007 


0.030 


0.005 


NGC6397 


0.024 


0.034 


0.007 


0.012 


0.133 


0.022 



Supplementary Table 3: Comparisons of the observed FWHMs of RGB stars and the 
measurement errors at the magnitude level of horizontal branch stars, V H b + 0.5 > V > 
V HB - 0.5. The measurement errors are those for individual stars and, therefore, mea- 
surement errors for individual subpopulation in GCs will be given by « a*(b - y)/y/n^ or 
a*(hk)/^/n^, where n p (> 20) is the number of stars in each subpopulation in GCs. Note 
that, while the FWHMs of most GCs have much larger values [> Qa*(hk)] than the mea- 
surement errors for individual stars in the hk index, the FWHM of NGC6397 is comparable 
in size to the measurement error in the hk index, consistent with the idea that NGC6397 
is the only normal GC in our sample (see Supplementary Figures I6l-f1~3~]). The last two 
columns, E(b - y)y 2 and E(hk) 1/2 , denote contributions due to the differential reddening 
effect assuming a 50% variation in the total interstellar reddening of each GC, with which 
observed FWHMs of GCs in the (b - y) color and in the hk index can not be explained 
simultaneously. 
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Supplementary Figure 1: (a) A comparison of Ca filter transmission functions between that in 
Anthony-Twarog et al. 12 (the black line) and that for the CTIO 1-m telescope (the blue line). Both 
filters have similar FWHMs, approximately 90 A, but the passband for the CTIO 1-m telescope 
is shifted approximately 15 A to the longer wavelength, (b) Synthetic spectra for the CN normal 
(the blue line) and the CN strong (the red line) RGB stars. The CN band at A 3885 A lies on the 
lower tail of the Ca filter but the contamination from the CN band is insignificant, (c) Synthetic 
spectra for the interstellar Ca II H & K lines. We adopt equivalent widths of 350 mA and 650 
mA for the interstellar Ca II H & K lines, respectively, with a gaussian line profile with a FWHM 
of 1 A (equivalent to Av r ~ 76 km/s). In the inset of the figure, the red line denotes the velocity 
profile for the interstellar Ca II H line and the blue line for the interstellar Ca II K line. This large 
amount of discrete interstellar absorption contributes only 0.010 mag to our results. 
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Supplementary Figure 2: (a & b) Blue crosses and red crosses denote RGB stars in the Ca-w and 
the Ca-s groups, respectively, with proper motion membership probabilities P > 90%. (c) RGB 
stars in the Ca-s group are shifted by Ahk = —0.20 mag to match with those in the Ca-w group, 
assuming the RGB split in M22 is due to differential reddening. The reddening correction value 
of Ahk = -0.20 mag for the Ca-s group is equivalent to E(B - V) = -1.69. (d) The RGB stars 
in the Ca-s group are shifted by A(b — y) = —1.25 mag, assuming E(b — y)/E(B — V) = 0.74, 
and two RGB sequences do not agree, in the sense that RGB stars in the Ca-s group is too hot to 
be in the RGB phase, (e & f) After applying reddening correction in V (= —5.24 mag), assuming 
Ay = 3. 1 x E(B — V) . The RGB stars in the Ca-s group become too bright to be members of M22, 
inconsistent with the proper motion study of the cluster. 
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Supplementary Figure 3: Spatial distribution of the Ca-w (blue dots) and the Ca-s (red dots) 
groups in M22. Note the absence of spatially patched features, indicating that differential redden- 
ing is not responsible for the RGB split in M22. 
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Supplementary Figure 4: Color-magnitude diagrams for M22 and two bulge fields (NGC6528 
and OGLEII-12). The black dots represent M22, the red dots and the blue dots denote NGC6528 
and OGLEII-12, respectively. The stars in the Milky Way bulge are fainter and redder than those 
in M22 are and they do not affect the double RGB sequences in M22. Black arrows indicate 
reddening vectors, assuming E(B — V) = 0.34 for M22 29 . 
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Supplementary Figure 5: (a) The color-magnitude diagram by Monaco et al. (b) Model 
isochrones for [Fe/H] = —1.6 (a blue line), and —1.4 (a red line), Y = 0.23, and 11 Gyr. (c) 
Model isochrones for [Fe/H] = -1.6, Y=0.23, 1 1 Gyr (a blue line) and [Fe/H] = -1.4, F=0.28, 1 1 
Gyr (a red line), (d) Model isochrones for [Fe/H] = -1.6, F=0.23, 1 1 Gyr (a blue line) and [Fe/H] 
= - 1 .4, F=0.28, 10 Gyr (a red line). 
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Supplementary Figure 6: The (b — y) and hk distributions of RGB stars in cu Cen. The red lines 
denote fiducial sequences of the cluster, which are forced to pass through the peak {b — y) colors or 
hk indices of given magnitude bins. The A(6 — y) and the Ahk are differences in the {b — y) color 
and the hk index, respectively, of each RGB stars from the fiducial sequences. The blue horizontal 
bars indicate measurement errors with a 2a* range (± la*) of individual stars at given magnitude 
bins. From the hk distribution, at least five distinct populations, whose hk splits are much larger 
than the measurement errors, can be found in u Cen. 
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Supplementary Figure 7: The (b — y) and hk distributions of RGB stars in M22. The blue 
horizontal bars indicate measurement errors with a 2a* range (± la*) of individual stars at given 
magnitude bins. Two distinct and discrete populations can be found in M22. At the magnitude 
of HB, the hk split between two populations is larger than 25 xa*(hk) or 250 x a p (hk), where 
a*(hk) and a p (hk) denote measurement errors for individual stars and populations in the hk index, 
respectively. 
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Supplementary Figure 8: The {b — y) and hk distributions of RGB stars in NGC1851. The blue 
horizontal bars indicate measurement errors with a 2a* range (± ler*) of individual stars at given 
magnitude bins. Two discrete populations can be found in NGC1851. At the magnitude of HB, the 
hk split between two populations is larger than 1 1 xa*(hk) or 55 xa p (hk). 
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Supplementary Figure 9: The {b — y) and hk distributions of RGB stars in NGC2808. The blue 
horizontal bars indicate measurement errors with a 2a* range (± la*) of individual stars at given 
magnitude bins. At least two discrete populations can be found in NGC2808. At the magnitude of 
HB, the hk split between two major populations is larger than 5xa*(hk) or 50 xa p (hk). 
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Supplementary Figure 10: The (b — y) and hk distributions of RGB stars in M4. The blue 
horizontal bars indicate measurement errors with a 2a* range (± ler*) of individual stars at given 
magnitude bins. Two discrete populations can be found in M4. At the magnitude of HB, the hk 
split between two populations is larger than 10 x a* (hk) or A5xa p (hk). 
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Supplementary Figure 11: The {b — y) and hk distributions of RGB stars in M5. The blue 
horizontal bars indicate measurement errors with a 2cr* range (± ler*) of individual stars at given 
magnitude bins. The RGB sequence of the cluster shows a large spread in the hk index, indicative 
of heterogeneous calcium abundances. At the magnitude of HB, the FWHM of RGB stars in M5 
is larger than 8 x a^hk) . 
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Supplementary Figure 12: The {b — y) and hk distributions of RGB stars in NGC6752. The blue 
horizontal bars indicate measurement errors with a 2<r* range (± lcr*) of individual stars at given 
magnitude bins. Two discrete populations can be found in NGC6752. At the magnitude of HB, the 
hk split between two populations is larger than \0xa*(hk) or 70xa p (hk). 
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Supplementary Figure 13: The {b — y) and hk distributions of RGB stars in NGC6397. The blue 
horizontal bars indicate measurement errors with a 2cr* range (± lcr*) of individual stars at given 
magnitude bins. It is the only normal GC in Figure 2. Note the similar degree of the RGB widths 
in A(b — y) and Ahk and the similar degree of RGB widths as the measurement errors. 



36 



0£) 

B 

m 
i 

I 



-3 
-2 

-1 h 

- 



1.0 - 



a 

<u 
a 



Z^ 0.5 

z 

0.0 



. 1 1 1 1 1 1 1 1 1 1 1 . 
■_ a j N2808 • 

++ + l .v + + + 

i | 1 >vH^pi"»r i | i 


i | i ill | i 

- f \ U4 - 

+ 

: M . .• : 
i i J 
+ +• + i * ++ + 
+ "i+ •'+-* 


1 | 1 1 II | 1 

- k ] .N675S ■ 
: : X "■ 

: ' "(!\- " 
+ .• .-r*V' 

: — , — | — | =i^^^ f " 


: b tf\ i 

. r1~ 1 1 ■ . 1 . " 


1 1 1,1 1 

\ JV 

^_ r+n *— 1 ,T , 1 


1 


* 1 i f] " 
l 

r ij 1 ~ 



-0.2 0.0 0.2 

Ahk (mag) 



-0.2 



2 
Z 



z 
z 



1.0 
0.5 
0.0 

1.0 
0.5 
0.0 



t — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 

- c 



0.0 

tshk (mag) 



-0.2 



o.o 



] — i — i — i — i — | — i — i — i — i — | — i — i — i — r 

^ h 



hhk (mag) 




-0.5 0.0 

[0/Fe] 



0.5 



-0.5 



0.0 0.5 
[O/Fe] 



-0.5 



0.0 0.5 
[O/Fe] 



- d 



- 1 — i — i — i — I — i — i — i — i — I — i — i — i — r - 




- - l 



_i i i 




- 1 — i — i — i — I — i — i — i — i — I — i — i — i — r~ 



- - n ,-, 



- 1 — i — i — i — i — i — i — i — r - 



-0.5 0.0 0.5 

[Na/Fe] 



-0.5 0.0 0.5 

[Na/Fe] 



-0.5 0.0 0.5 

[Na/Fe] 




l 1 1 1 1 1 1 r 



J 



- - o 



r n H 

111 ' I 




"T" 



- 1 1 r- 



_1_ 



[Fe/H] 



-1.2 -1.0 

[Fe/H] 



-1.6 -1.4 

[Fe/H] 



Supplementary Figure 14: (a) A plot of V - V HB versus Ahk for RGB stars in NGC2808. The 
blue and red plus signs denote the Ca-w and the Ca-s RGB stars. The dashed line denotes the 
boundary between the two groups at Ahk = —0.05 mag. (b) The Ahk distribution of NGC2808 
RGB stars, (c) The [O/Fe] distributions of the two RGB populations in NGC2808. The shaded 
histogram outlined with blue color is for the Ca-w group and the blank histrogram outlined with 
red color is for the Ca-s group, (d) The [Na/Fe] distributions, (e) The [Fe/H] distributions, (f-j) 
Same as (a-e), but for M4 RGB stars with the boundary at Ahk = —0.08 mag. (k-o) Same as 
(a-e), but for NGC6752 RGB stars with the boundary at Ahk = -0.04 mag. 
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Supplementary Figure 15: (a) A plot of V — Vhb versus Ahk for 40 RGB stars in u Cen of 
Johnson et al. 5 (b - e) Elemental abundances of 40 RGB stars in uj Cen as functions of Ahk. (f - g) 
Comparisons of our photometric metallicity, [Fe/H]^, with spectroscopic metallicity, [Fe/H] spec , 
and calcium abundance, [Ca/H] spcc . The linear fits to the data are shown with red lines. 
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Supplementary Figure 16: Metallicity distribution functions for eight GCs derived from the hk 
index. In the figure, [Fe/H] hk,corr is our recalibrated photometric metallicity using the equation 
d3]). Note that we only use bright RGB stars in order to minimize contamination from off-cluster 
field and red-clump populations. For most GCs, signs of multiple stellar populations persist in our 
MDFs. 
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Supplementary Figure 17: Comparisons of synthetic spectra for [Fe/H] = —1.6, Teff = 4750 K, 
log g = 2.0. (Upper panel) The red line denotes synthetic spectrum for [Ca/Fe] = 0.25 dex and the 
blue line denotes that for [Ca/Fe] = 0.45 dex. We adopt the fixed aluminium abundance of [Al/Fe] 
= 0.50 dex. (Lower panel) The red line denotes the synthetic spectrum for [Al/Fe] = 0.00 dex and 
the blue line denotes that for [Al/Fe] = 1.00 dex. We adopt the fixed calcium abundance of [Ca/Fe] 
= 0.30 dex. The effect of aluminium contamination on the hk index appears to be negligible. 
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